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Using a new instrument, the HEAT collaboration has confirmed the excess of cosmic ray positrons 
that they first detected in 1994. We explore the possibility that this excess is due to the annihilation 
of neutralino dark matter in the galactic halo. We confirm that neutralino annihilation can produce 
enough positrons to make up the measured excess only if there is an additional enhancement to 
the signal. We quantify the 'boost factor' that is required in the signal for various models in the 
Minimal Supersymmetric Standard Model parameter space, and study the dependence on various 
parameters. We find models with a boost factor > 30. Such an enhancement in the signal could arise 
if we live in a clumpy halo. We discuss what part of supersymmetric parameter space is favored (in 
that it gives the largest positron signal), and the consequences for other direct and indirect searches 
of supersymmetric dark matter. 



I. INTRODUCTION 

Several years ago the HEAT collaboration reported an 
excess of cosmic ray positrons with energies ~ 10 GeV 
0. In the past year they again measured this excess 
using a new instrument, and found excellent agreement 
In this paper we revisit the possibility that this ex- 
cess is due to annihilations of Weakly Interacting Massive 
Particles (WIMPs) in the galactic halo, in particular the 
neutralinos in supersymmetric models. The possibility of 
detecting positrons as annihilation products of WIMPs 
in the galactic halo has been discussed previously, both 
as a continuum and as a monochromatic source |^-||. 
We update the results of including the new data, and 
confirm the conclusions of that paper, finding that neu- 
tralino annihilation can not produce enough positrons to 
make up the measured excess without an additional en- 
hancement to the signal. Recently, this point has been 
reiterated in Ref. 

In calculating the observed positron flux from annihila- 
tions in the halo, we encounter several astrophysical un- 
certainties. First, cosmic ray propagation is not perfectly 



understood, though the errors are unlikely to be larger 
than a factor of two. More importantly, the structure of 
the galactic dark halo is unknown. Any dumpiness in the 
halo serves to enhance the signal, whether it is a single 
nearby clump (or one containing the Earth) , or a uniform 
distribution of clumps. There is no compelling argument 
for any particular value of the enhancement factor, be 
it unity or in the thousands or more. In this paper we 
carefully discuss the possibility that a clumped galactic 
halo could account for the measured positron excess. 

We define Bs to be the boost factor that the WIMP 
annihilation signal from a smooth galactic halo must be 
multiplied by to match the HEAT data. We have ex- 
plored models in the Minimal Supersymmetric Standard 
Model (MSSM) parameter space to find how large the 
boost factor must be for each of the models. The low- 
est boost factor we found is roughly 30 for a WIMP 
that is primarily a Bino in content with mass 160 GeV. 
For Bs < 100, we find that the models are gaugino- 
dominated, though some have significant Higgsino frac- 
tions. The masses of the models are in the range 150-400 
GeV for the most part. For 100 < Bg < 1000, the masses 
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are as large as 2 TeV, and some very pure Higgsinos be- 
come allowed. For both cases there are a significant num- 
ber of models which have a large contribution to so 
that the anomalous magnetic moment of the muon can 
be explained. We have investigated the dependence of 
the boost factor on various parameters. There is essen- 
tially no dependence on parameters such as tan (3 or toq 
(the sfermion mass scale) . The boost factor does depend 
strongly on the relic density of WIMPs; the lowest boost 
factors are required for the models with the smallest relic 
density without rescaling (defined below). 

One can ask the question: even with a large boost 
factor, can neutralino annihilation produce the "bump" 
seen in the positron spectrum just above 10 GeV? Even 
if there is a line signal direct from the annihilation, it 
gets spread out by the propagation, so that a bump does 
not get produced. Thus, if one looks "by eye" , one would 
conclude that neutralinos cannot produce the data. How- 
ever, this is an inappropriate way to ask the question. 
One should instead study the problem statistically to see 
if one can find a neutralino model with a good • With 
a combination of background and annihilation signal, wc 
are able to find statistically reasonable fits to the spectral 
shape for boost factors above 30. 

II. SUPERSYMMETRIC MODEL 

We work in the Minimal Supersymmetric Standard 
Model (MSSM). In general, the MSSM has many free 
parameters, but with some reasonable assumptions we 
can reduce the number of parameters to the Higgsino 
mass parameter /i, the gaugino mass parameter M2, the 
ratio of the Higgs vacuum expectation values tan /?, the 
mass of the CP-odd Higgs boson tua, the scalar mass pa- 
rameter mo and the trilinear soft SUSY-breaking param- 
eters Ai, and At for the third generation. In particular, 
we don't impose any restrictions from supergravity other 
than gaugino mass unification. For a more detailed defi- 
nition of the parameters and a full set of Feynman rules, 
see Refs. [|o[||. 

The lightest stable supersymmetric particle in the Min- 
imal Supersymmetric Standard Model (MSSM) is most 
often the lightest of the neutralinos, which are superposi- 
tions of the superpartners of the neutral gauge and Higgs 
bosons, 

xl^NiiB + Ni2W'^ + Ni^H^ + NuH^. (1) 

For many values of the MSSM parameter space, the 
relic density Q.^h? of the (lightest) neutralino is of the 
right order of magnitude for the neutralino to consti- 
tute at least a part, if not all, of the dark matter in 
the Universe (for a review see Ref. Q). Here VL^ is 
the density in units of the critical density and h is the 
present Hubble constant in units of 100 km s"^ Mpc^^. 
Present observations favor h = 0.7 ± 0.1, and a total 
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-3 -3 
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3 3 



TABLE I. The ranges of parameter values used in the 
MSSM scans of Refs. |io[0,^jl^,§|l|,0 . 



matter density Qm = 0.3 ± 0.1, of which baryons con- 
tribute roughly Qbh'^ « 0.02 @. Thus we take the 
range 0.05 < O^/i^ < 0.25 as the cosmologically inter- 
esting region. This region can be narrowed somewhat if 
we consider the results of CMB anisotropy measurements 
(summarized in e.g. |jl^), which favor Q^h^ = 0.14±0.05. 
We are also interested in models where neutralinos are 
not the only component of dark matter, so we separately 
consider models with arbitrarily small fl^h^. 

As a scan in MSSM parameter space, we have used 
the database of 

MSSM models buih in Refs. 0,|ll],^6|jl|,|,|l|,|l9| . The 
overall ranges of the seven MSSM parameters are given 
in Table ||. While the ranges are extreme, most in- 
teresting models fall in a much more modest region of 
parameter space, with the notable exception that very 
pure Higgsinos, thus very large M2, can not be ruled 
out at present. The database embodies one-loop cor- 
rections for the neutralino and chargino masses as given 
in Ref. ||2^, and leading log two-loop radiative correc- 
tions for the Higgs boson masses as given in Ref. 111. 
For all of the MSSM models in the scan of parameter 
space, the database contains results for expected detec- 
tion rates of the particles in a variety of neutralino dark 
matter searches. The database includes the relic den- 
sity of neutralinos fi^/i^. The relic density calculation 
in the database is based on Refs. (l^,^^ and includes 
resonant annihilations, threshold effects, finite widths of 
unstable particles, all two-body tree-level annihilation 
channels of neutralinos, and coannihilation processes be- 
tween all neutralinos and charginos. The database also 
includes the supersymmetric correction to the anomalous 
magnetic moment of the muon = {g^ — 2)/2 which 
is important for dark matter searches in light of new 
data | p3[ | indicating a deviation from the standard model 
prediction, as discussed by e.g. Ref. |2j]. In this paper 
we will identify models which have a large contribution, 
10 X 10-1° < Aa^(SUSY) < 75 x 10"!°, to the anoma- 
lous magnetic moment of the muon as being particularly 
interesting. 

We examined each model in the database to see if it is 
excluded by the most recent accelerator constraints. The 
most important of these are the LEP bounds HI on the 
lightest chargino mass 

r 88.4 GeV , |m^+ - m^o| > 3 GeV 

\ 67.7 GeV , otherwise, ^ ' 

and on the lightest Higgs boson mass rrih (which ranges 
from 91.5-112 GeV depending on tan/3) and the con- 
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straints from b — > 57 
in DarkSUSY [|7|). 



(we use the LO implementation 



III. POSITRON FLUX 

We obtain the positron flux from neutraHno annihila- 
tion in the galactic halo following Ref. 0. The model 
is a true diffusion model and assumes that the diffu- 
sion region of tangled galactic magnetic field is an infi- 
nite slab. This approximation is reasonable since most of 
the positrons are emitted quite nearby so that the outer 
radial boundary is unimportant. Furthermore, energy 
losses due to synchrotron radiation and inverse Comp- 
ton scattering from the cosmic microwave background 
and from starlight are included. This model roughly 
agrees with earlier work [Q, though the inclusion of in- 
verse Compton scattering from starlight is crucial as it 
doubles the energy loss rate. 

As we will discuss in the following sections, the 
positron fiux from a smooth galactic halo is too low to 
explain the positron excess, as has been discussed previ- 
ously However, any deviations from smoothness 
serve to enhance the annihilation signal, as the anni- 
hilation rate is proportional to the neutralino density 
squared. However, we must be careful that in postu- 
lating a boost factor, we do not overproduce the other 
products of neutralino annihilation, especially antipro- 
tons and gamma rays ||2^. We do have some freedom 
here, in that the boost factors for positrons, antipro- 
tons and gamma rays are not necessarily equal, as their 
propagation is not the same. For example, a nearby 
clump would serve to increase the positron fiux more than 
it would increase the antiproton or gamma ray fiuxes, 
as positrons have the shortest range (they have shorter 
range than the antiprotons because of their rapid energy 
loss) . Many of the antiprotons come from far away, out- 
side the location of the clump, almost a third of them 
from as far away as the Galactic Center. Positrons, on 
the other hand, come from much closer, roughly within 
a few kiloparsecs. 

We will fit the full positron dataset of the HEAT ex- 
periments (1994 and 1995 combined data ||l[ and the 
2000 data j2|). We use the positron fraction data, as 
the error bars are smaller and the data cleaner. The full 
dataset consists of twelve independent measurements of 
the positron fraction at various energies. 

We will in the following assume that the standard pre- 
diction for the positron background |2^,^ is correct to 
within a normalization factor A''. We are aware that cos- 
mic ray propagation is not completely understood, and 
that even the best efforts to reproduce the observed cos- 
mic ray spectra need to rely on yet-to-be-understood ad 
hoc assumptions on the dependence of the diffusion con- 
stant on energy and on the source spectrum pO| ]. How- 
ever, we gather from the latter work that the discrepan- 
cies between the observed and theoretical positron spec- 



tra lie preponderantly at energies smaller than a few 
GeV, where they can become as large as a factor of 4 
in the hundreds of MeV range. At the slightly higher en- 
ergies where the HEAT bump is, the theoretical models 
in [ ^0| , which cover a wide range of theoretical assump- 
tions, agree to within 20%. While this may give some 
justification to our use of model 08-005 of Ref. ||2^ as 
our standard positron background, we nevertheless stress 
that it may be possible to explain the positron bump by 
purely astrophysical means (although we do not know 
how). Keeping these uncertainties in mind, we proceed 
with the assumption that the background calculation is 
correct, and we study the possibility that neutralino an- 
nihilation can account for the excess positrons. 

We assume that the positron signal from neutralino 
annihilation can be rescaled by a normalization factor 
(boost factor) Bs- We find that the best fit normaliza- 
tion of the background with no signal from neutralinos is 
N = 1.14, with — 3.33 per degree of freedom. When 
adding the signal, we make a simultaneous fit of the nor- 
malization of the background N and the normalization 
of the signal Bs, for each supersymmetric model in the 
database. We say that a given model "gives a good fit to 
the positron data" when: (1) the background-plus-signal 
fit fits the data better than the background-only fit with 
a decrease in per degree of freedom greater than unity, 
namely the background-plus-signal fit has ^ 2.33 per 
degree of freedom; (2) the best fit normalization of the 
background N is between 0.5 and 2.0, namely the calcu- 
lation of is correct to within a factor of two according 
to the best fit. 

The positron fiuxes are more than an order of mag- 
nitude smaller than the HEAT measurements, and we 
find that the best fit normalizations of the signal Bg lie 
between 30 and 10^°. Values of Bg as large as 10^" are 
hardly realistic, but Bs up to 100-1000 might be accept- 
able given the uncertainties in the halo structure (the 
halo could e.g. be clumpy 28 ,3^]). 



A. Antiprotons 

In addition, we require that the antiproton fiux from 
annihilations (3^jl^ not be too large, given the boost 
factor required for each model. There is a significant 
correlation between the antiproton and positron fiuxes 
due to neutralino annihilations (see e.g. Fig. 8 of Ref. [Q), 
so this constraint is important. Following Ref. we 
take the antiproton fiux as 



$p = fc(l + 0.75Bs)^' 



smooth 



(3) 



where k represents the difference in enhancement fac- 
tors between the antiprotons and positrons. The factor 
0.75 comes from the fact that the antiprotons that reach 
the Earth on average are produced further away than 
the positrons. In particular, the antiprotons produced 
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close to the galactic center make a significant contribu- 
tion to the flux at Earth. In these denser environments, 
a clumpy distribution enhances the signal less than in 
our local environment, and hence the different scalings 
of the positron and antiproton fluxes (see Ref. for 
more details). 

We take the constraint from the combined 1995 and 
1997 data of the BESS collaboration Q. Note that we 
use the BESS data rather than the HEAT 2000 antipro- 
ton data (the new instrument measures both positrons 
and antiprotons) for the following reason: the observed 
antiproton flux rises with momentum to a maximum 
value and then falls. The spectrum from neutralinos, on 
the other hand, is flat with momentum and then cuts off 
before it reaches the observed peak. Hence the strongest 
constraint on neutralinos comes from lower energy mea- 
surements. BESS goes to much lower energy than HEAT 
2000 and thus places a stronger constraint. The BESS 
collaboration measured the cosmic ray antiproton flux at 
low energies to be 

$p(r = 400 - 560 MeV) = 1.27^1^^2 

X 10"^(cm2 s sr GeV)-^ (4) 

Using the central value as the maximum allowed anni- 
hilation flux, and taking k — 1, we find no models with 
a boost factor Bg < 100, though there are a handful 
of models with Bs < 300, including several with a sig- 
nificant value of a^. Taking k — 0.2, we find the con- 
straint much less punishing, and given the uncertainties, 
we choose this value instead. 



B. Two Successful Models 

In Fig. |l| we plot the positron data from the HEAT 
94-1-95 and HEAT 2000 experiments, together with the 
background only fit, and two interesting SUSY models 
that have good fits as well as large contributions to a^. 
The antiproton constraints have been applied with k — 
0.2. Note that we have found the boost factor to be only 
weakly dependent on tan/3, so that the values of tan/9 
in the figure do not play any important role (except that 
tan/3 is correlated with the SUSY contribution to a^). 
For other models see Fig. 7 of Ref. §. 

The apparent sharp increase in the positron fraction 
around 7 GeV is not evident in any of our SUSY mod- 
els, even before the smoothing effects of energy loss on 
the spectrum. In principle, positrons from direct gauge 
boson decays have a perfectly flat spectrum (before prop- 
agation) with cutoffs at 




where m is the gauge boson mass. For W , a fea- 
ture at 7 GeV would thus be had for ~ 238 GeV. 



However there are also positrons from hadronizations at 
least from the hadronic gauge boson decays, and possibly 
from direct annihilations to quark-antiquark pairs. The 
hadronic component is dominant at the lower cutoff (but 
not always at the upper cutoff), which means that we 
can not reproduce a sharp bump at 7-8 GeV as indicated 
by the data, but rather a smoother bump over a larger 
energy range as seen in Fig. |l|. 

A way to sharpen the neutralino annihilation positrons 
into a bump is to have them all come from a nearby clump 
which is smaller than the propagation length. Then a line 
signal would not be smeared out. This problem has not 
yet been treated in depth. It requires a different solution 
of the diffusion equation and is the subject of a future 
work. 



IV. FAVORED REGION IN SUPERSYMMETRIC 
PARAMETER SPACE 

Having computed the positron flux and required en- 
hancement factors to give a good fit to the positron data, 
we can now study the supersymmetric parameter space 
and identify the favored regions. The composition of the 
neutralino, namely if it is gaugino or Higgsino, is per- 
haps its most interesting property. As our indicator of 
composition, we use the gaugino to Higgsino ratio 

Zg ^ |jVll|' + |iVl2|' , . 

1-Zg |/Vi3P + |/Vi4|2- ^' 

In Fig. ^ we plot this ratio vs. the neutralino mass sepa- 
rately for models with Bs < 100 and 100 < B^ < 1000, 
good J good background normalization as discussed in 
the previous section, and also a relic density in the region 
favored by the CMB, fl^h^ = 0.14 ±0.05. For B^ < 100, 
we find that the models are gaugino-dominatcd, though 
some have significant Higgsino fractions. The masses of 
the models are in the range 150-400 GeV for the most 
part. For 100 < < 1000, the masses are as large as 
2 TeV, and some very pure Higgsinos become allowed. 
For both cases there are a significant number of models 
which have a large contribution to a^. 

V. BOOST FACTOR 

It is instructive to study the best fit boost factor as a 
function of the supersymmetric parameters of the models 
under discussion. As in the previous sections, we will 
restrict ourselves to models that provide a good fit to the 
HEAT data, have a relic density in line with CMB data, 
and do not produce an overabundance of antiprotons. 

We first study the dependence of the relic density on 
the boost factor in these models. To do so we will of 
course neglect to apply the constraint on relic density, 
though we retain all other constraints. We will show 
two cases, first assuming that the dark halo density is 
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FIG. 1. Positron fraction data and fits. We illustrate positron data from HEAT 94+95 and HEAT 2000, a background only 
fit, and a SUSY+background fit from two interesting models from the MSSM database. Two additional curves separately 
display the SUSY and background components of the combined SUSY+background fit. These models are gaugino dominated 
and have contributions to in line with the experimental discrepancy. The model in Fig. la has positrons primarily from 
hadronization, while the model in Fig. lb has hard positrons from direct gauge boson decays. 
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FIG. 2. Neutralino composition vs. mass for well-fitting models. For Bs < 100 these are mixed and gauginos, mostly from 
150 — 400 GeV. For 100 < Bs < 1000 the masses extend to 2 TeV, and very pure Higgsinos are also allowed. In both cases 
many of the models have a contribution to in line with the measured discrepancy. 
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independent of the relic density, and second applying a 
rescaling. This rescaling is applied for models whose relic 
density is less than fl^h"^ = 0.09, the lower value in the 
CMB range, and is defined as follows. For low relic den- 
sities, neutralinos would only make a fraction of the dark 
halos of galaxies, and in principle that fraction should be 
proportional to the relic density. 



CDM 



PCDM,gal. 



(7) 



Here, subscript "gal" indicates that the density is that 
inside the Galaxy and subscript "CDM" refers to the 
dominant matter component. As annihilation depends 
on the square of the density, we rescale as the square 
of the fraction. The rescaling affects the best fit boost 
factor and the antiproton flux as follows: 



$55 

p 



B, 



smooth 



0.09 



^smooth / ^X^^ 
0.09 



(8) 
(9) 



In Fig. 1^ we plot the boost factor versus the relic den- 
sity for both cases, not rescaled and rescaled. Without 
rescaling, there is a clear trend that Bg increases linearly 
with ^^h? . This is expected because the relic density is 
inversely proportional to the annihilation cross section, 
and so is the boost factor. When taking the rescaling 
into account, we find that the lowest boost factors are 
required for the models with the smallest relic density 
without rescaling, that is with Vlh^ = 0.09 according to 
our choice. 

We take two main points from Fig. ^. First, our pre- 
ferred region depends somewhat sensitively on the cuts 
we make in the relic density. Enlarging our definition 
of the region of cosmological interest would have a sig- 
nificant effect on the number of allowed models in the 
database. Our conclusions would however be broadly 
similar. Second, and more importantly, we see the fun- 
damental problem of explaining the positron excess with 
neutralino annihilation. Given the observed value of 
the dark matter density, the expected annihilation cross 
section is too small to explain the observed excess of 
positrons without some boost due to clumping or some 
other mechanism (for example some models discussed in 
Refs. 35 1). This point is independent of the specific 
model for the WIMP, and only relics on the fact that the 
relic density is due to a thermal freeze-out of a stable 
(or very long-lived) species, and reasonable annihilation 
branching fraction to hadrons. 

Concerning the other supersymmetric parameters, we 
find that the boost factor is only weakly dependent on 
tan/3 and mo, the sfermion mass scale. We see a rough 
trend that models with heavy neutralinos need larger 
boost factors, but this is simply related to the fact that 
the number density scales as the inverse of the mass and 
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FIG. 4. Boost factor versus antiproton flux. The trend that 
models with small antiproton fluxes require large boost factors 
in the positron signal reiterates the statement that there is a 
signiflcant correlation between the antiproton and positron 
fluxes. The hatched band indicates the BESS measurement 



[0- 



thus the annihilation rate scales as the inverse square of 
the mass. 

Since the constraint from the antiproton flux is so im- 
portant, we now show how this quantity depends on the 
required boost factor. Of course we now neglect to apply 
the constraint on antiproton flux, but we retain all other 
constraints. In Fig. ^ we plot the boost factor vs. the an- 
tiproton flux in the 400-560 MeV bin for easy comparison 
with the BESS experiment, shown as the hatched band. 
That small antiproton fluxes imply large boost factors 
is another statement of the fact that the antiproton and 
positron fluxes are significantly correlated. Furthermore, 
we see the advantage of allowing fc < 1, as the bound on 
antiproton flux sits at 1.27loj2 x 10"^ (cm^ s sr GeV)"^ 
Even for k — 0.5, a significant number of models becomes 
allowed, especially with boost factors Bg < 300. 

Finally, we comment on the feasibility of the required 
boost factors, Bg > 30- It is well known that dark matter 
is clumpy in a large range of length scales; such clumps 
are clusters, galaxies, dwarf galaxies, etc. On such large 
scales the enhancement factor is in excess of 100 accord- 
ing to simulations of large scale structure |Q . The ques- 
tion for us is whether such clumps persist at scales smaller 
than several kiloparsecs, which is the size of the emission 
region for positrons detectable at the Earth. Unfortu- 
nately, there is really no data at these distance scales, 
either observational or from simulations. Without evi- 
dence to the contrary, we must allow such enhancements 
to be possible. 
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FIG. 3. Boost factor versus relic density. In the first panel, no rescaling is done. The trend that small boost factor indicates 
small relic density is clearly seen. In the second panel, the rescaling is performed for models whose relic density is less than 
Q,-^h^ — 0.09, and it is clear that the smallest boost factors come from the smallest relic densities that require no rescaling. 



In obtaining a boost factor, we have assumed that we 
can average over a volume containing many small clumps. 
If the halo is not smooth, but we can average over a large 
volume (relative to the propagation length of positrons 
from several clumps), then we can pull out an enhance- 
ment factor. In other words, we can use the results of 
the DarkSUSY code for a smooth halo and just multiply 
the result by a boost factor. 

VI. OTHER DARK MATTER SEARCHES 

In order to be convinced of an exotic interpretation 
of cosmic ray data, we would like confirmation by some 
other technique. In this section we discuss other dark 
matter search techniques that might give us more confi- 
dence that the positron excess really is due to an exotic 
primary component. In particular, we discuss direct de- 
tection of neutralinos by elastic scattering, indirect de- 
tection by gamma ray lines, and furthermore by neutri- 
nos from capture and annihilations in the centers of the 
Earth and Sun. 



A. Direct detection 

Direct detection of galactic halo neutralinos is one of 
the most promising techniques for detecting dark mat- 
ter, and there are several experimental collaborations 
undertaking this program, e.g. DAMA CDMS 



CRESST ||], EDELWEISS lojLCryoarray M, GE- 
NIUS IGEX (CanFranc)|g2|, HDMS |3[, MI- 
BETA, ROSEBUD m, LiF/TOKYO ||, UKDMC 
!§, SACLAY, ELEGANT V |47|, and Baksan. In the 
next ten years it is expected that neutralinos with elastic 
scattering cross sections on nucleons as low as 10~^ pico- 
barns or perhaps even lower can be probed |^,^. The 
rates in detectors only depend on the local halo density 
at present, so they will not put any severe constraints on 
the dumpiness of the halo as a whole. These rates can 
of course be enhanced if we happen to be inside a clump 
at present [Q. 

Direct detection is especially exciting in light of the 
measurement of a possible discrepancy in the anomalous 
magnetic moment of the muon |23|. Models with large 
contributions to tend to also have large elastic scat- 
tering cross sections [ p^ , and we find that many of the 
models that can explain the positron excess with boost 
factors Bs < 1000 also have large contributions to a^. 
We plot the scattering cross section in Fig. |[ 

B. Neutrinos from the Earth and Sun 

Another possible method to detect neutralino dark 
matter is neutrino telescopes, such as at Lake Baikal [ p9[ , 
Super-Kamiokande |5^], in the Mediterranean |5|], and 
the South Pole Neutralinos in the galactic halo un- 
dergo scatterings into bound orbits around the Earth and 
Sun, and subsequently sink to the centers of these bodies, 
possibly giving a detectable annihilation signal in neutri- 
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FIG. 5. Neutralino-nucleon elastic scattering cross section. 
We only show models passing the goodness of fit cuts and 
require Bs < 1000. The current CDMS exclusion ^ (solid 
line) is plotted along with the expected reach of the CDMS, 
CRESST H (dotted line) and GENIUS (dashed line) 
experiments. 



nos at GeV and higher energies [Q. The detectabihty 
of this signal is strongly correlated with the neutralino- 
nucleon cross section. As discussed in Ref . , the signal 
is not likely to be sensitive to the dumpiness in the halo. 
This statement assumes that equilibration time between 
capture and annihilation is relatively long and that, av- 
eraging over the lifetime of the Galaxy, the average local 
density is neither overdense nor underdense. 



C. Gamma rays 

Gamma rays from annihilations, both a continuum 
component [p4[ and a monochromatic component fs^ ] 
may provide another handle on neutralinos in the galac- 
tic halo. Experiments such as the GLAST satellite js^ ] 
and Atmospheric Cerenkov Telescopes (ACTs), such as 
VERITA S i7| , STACEE |§, CANGAROO-III and 
MAGIC ||60|] , may have the necessary sensitivity to detect 
annihilation photons in our galaxy above the background 

0- . . . 

To minimize the impact of the halo model and of exper- 
imental uncertainties, we concentrate on the flux at high 
latitudes, b > 60° and 0° < £ < 360° (Afi = 0.84 sr), 
although we also consider the flux towards the galactic 
center. A modified isothermal profile gives 



J(90°) = 0.93(1 + 1.8Bs), 
where the gamma ray flux is given by 



(10) 



$^ 1.878 X 10"^^ (cm^ s sr 

N^{av) 



100 GeV 



X 
-2 



(11) 



There is only a very weak halo model dependence in this 
result for J at high galactic latitude We might ex- 
clude models which have too high a gamma ray Hux as 
compared with the measured value at high latitude 61|, 



$^(£; > 1 GeV) = (1.0 ± 0.2) x 10"^(cm2 s sr) 



(12) 



though with boost factors Bg < 100 the antiprotons are 
always more powerful p8| . However, boosting the signal 
of the gamma ray lines may allow their detection, which 
would be a clear confirmation of the neutralino halo. 

The sensitivities of gamma-ray detectors to the gamma 
ray lines can be computed following Ref. First the 
exposure is determined as a function of energy, as the 
ACTs in particular have an effective collection area that 
depends on energy. For the ACTs we consider, these 
are of order 10^ cm^ near threshold and rising to 10^ 
and more at TeV energies. ACT integration times of 500 
hours are assumed, while a 2 year GLAST integration is 
assumed. The GLAST exposure is taken to be a constant 
1800 cm^ sr, which simply multiplies the 2 yr integration, 
the fraction of time pointing towards a target already ac- 
counted for. The angular field of view for the ACTs is 
taken to be 0.01 sr, a circle 3.5° in radius. Based on 
these exposures, the number of background events is de- 
termined from the extragalactic gamma ray background, 
and additionally for ACTs, the backgrounds of cosmic 
ray electrons and misidentified hadrons. In fact the pho- 
ton background is unimportant for ACTs: 

$had = 1.0xlO-M__j (cm^ssrGeV)"', (13) 
= 6.9 X 10-^ f j (cm^ s sr GeV) " , (14) 

= 6.0 X 10"^ f j (cm^ ssrGeV)~\ (15) 

We note that the background flux for a gamma ray line at 
a specific energy need only be integrated over the energy 
resolutions of the experiments, taken to be fractionally 
0.15 for ACTs and 0.015 for GLAST. 

The signal is obtained from Eq. ([ll|), taking care to 
properly average J over the field of view. For high galac- 
tic latitudes and for ACTs towards the galactic center 
this is a minor consideration as J changes little over the 
field of view, though for GLAST we find that at the galac- 
tic center the averaged value is about ten times smaller 
than the central value. We require a 5cr excess above 
background to claim a detection. 

In Fig. 1^ we plot the flux from the gamma ray lines at 
high latitude for models with Bg < 1000, appropriately 
boosted by Bg. We include the expected sensitivity of 
the VERITAS l57tl and MAGIC M ACTs as well as the 
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neutralino mass (GeV/c^) 

FIG. 6. Gamma ray line flux. We include the expected 
sensitivityfsolid lines) of the VERITAS @, MAGIC @ and 
GLAST |56| experiments. The left axis is the flux from high 
galactic latitude (dotted sensitivity curves), and the right axis 
is the flux from the galactic center for an isothermal halo with 
a 1 kiloparsec core (solid sensitivity curves). We show only 
models passing the goodness of flt cuts and require Ba < 1000. 

GLAST satellite. Furthermore, we include a pre- 
diction appropriate to the galactic center, assuming an 
isothermal halo with a 1 kiloparsec core, which has a sig- 
nal roughly 1000 times larger, J ~ 1600i3s (this value 
is decreased by a factor of 10 when averaged over the 
GLAST field of view). At high latitude, a more sen- 
sitive experiment is probably required, though towards 
the galactic center, many models would give detectable 
fluxes in the gamma ray lines. 

VII. CONCLUSIONS 

The cosmic ray positron excess is intriguing, as there 
is no simple astrophysical model that can explain it. We 
are left to consider a primary component, such as from 
neutralino annihilations. We summarize here our conclu- 
sions concerning this scenario. 

First, the observed value of the dark matter density 
implies (assuming thermal production) an annihilation 
cross section that is too small to reproduce the positron 
excess without some form of enhancement. This is a gen- 
eral statement, not tied to a specific model. We thus 
resort to enhancing the signal; fortunately such an en- 
hancement is natural as the dark halo is expected to be 
clumpy. This leads to a second difficulty, namely that 
one can not enhance the positron signal without enhanc- 
ing other signals, especially antiprotons which would also 



be produced by annihilation. Hence antiprotons provide 
a further constraint on this scenario. Indeed lowering the 
antiproton flux is a further, albeit small, price to pay in 
the neutralino annihilation scenario. (Note that antipro- 
tons come from much farther away than the positrons, so 
their fluxes are not always directly correlated.) In addi- 
tion, in order to obtain a positron spectrum that matches 
all the data, we had to adjust the normalization of the 
background (another price we had to pay). We had to 
choose a positron background a factor of 2 lower than 
the standard fit to the positron data with background 
alone. The reason for this lowered background normal- 
ization is that one cannot overshoot the data at energies 
1-3 GeV. In reality the background is not terribly well 
understood, and though it cannot by itself explain the 
turnup in the data at 10 GeV, one wonders if perhaps the 
boost factor might not be plausibly lower than the values 
we have found. However, we find this possibility unlikely. 
We should mention here that the propagation uncertain- 
ties make the change in background normalizations and 
relative boosts between antiprotons and positrons (the 
A: = 0.2vs. fc = l issue) more plausible, and we do not 
believe these to be serious concerns with our analysis. 

Second, assuming that the boost factor is between 30 
and 100, we find gaugino-dominated SUSY models that 
satisfy all constraints, have neutralino masses in the 150- 
400 GeV range, and have a large contribution to the 
anomalous muon magnetic moment a^. Allowing boost 
factors as large as 1000 extends the mass range to 2 TeV, 
and furthermore allows Higgsino-dominated neutralinos. 
Such boost factors are certainly plausible, and with no 
evidence to the contrary, we must take this possibility 
seriously. 

Confirmation of the annihilation hypothesis could 
come from several approaches. The direct detection of 
halo neutralinos would certainly be a powerful indica- 
tor, as would neutrinos from the center of the Earth and 
Sun. Antiprotons and gamma rays could help study the 
dumpiness of the galactic halo, helping to determine if 
it is in fact as large as in the scenario we have presented. 
In particular, boosting the intensity of gamma ray lines 
may allow their detection, which would be a clear confir- 
mation of a neutralino halo. The next several years will 
be an exciting time for particle dark matter searches. 
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